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SUMMARY 


12 


Proper use of reduced data from falling sphere soundings 
requires a detailed knowledge of the many factors affecting their 
accuracy. Criteria for judging the reliability of ROBIN sphere 
data (wind and density) are reviewed, and the effect of the data 
smoothing interval on the suppression of radar tracking error and 
the retention of real atmospheric detail is briefly considered. 

The height range of valid temperatures is described in relation to 
the prevailing pressure scale height and the climatological temper- 
ature regime. ROBIN soundings considered valid by stated criteria 
are presented which illustrate their fidelity to large-scale atmo- 
spheric variations . Also shown are rare examples of soundings 
satisfying the above criteria which indicated physically improbable 
variations. Special data comparisons include Arcasonde and ROBIN 
temperatures from Ascension Island, which indicate fair agreement 
near the stratopause; and densities based on thermistor and sphere 
data. The latter reveal a systematic difference of about 8%; var- 
ious explanations are considered, but none is found to account for 
the full difference. The article concludes with a brief review of 
past uses of sphere data, together with remarks on persisting pro- 
blems relevant to all sphere measurements. 


INTRODUCTION 


Perhaps no method for sounding the upper atmosphere with roc- 
kets has required so much study as the falling sphere technique. 

Much of the commentary has concerned the ROBIN sphere (Engler, 

> Luers , ref. 2) , but some of the problems of interest apply 
to other spheres as well (Jones and Peterson, ref. 3; Salah, ref. 4; 
Faucher et al., ref. 5; Peterson, ref. 6; Champion and Faire , ref. 7). 
Nearly 700 ROBIN soundings have been obtained, including a research 
series of 188 observations in 1960-62 (Lenhard and Kantor, ref. 9) 
and shorter research series in 1965 and 1966. More than 100 sound- 
ings have been taken with other spheres, and these have a propor- 
tionally greater value because of the higher altitudes reached. 
Densities from ROBIN soundings have been obtained generally in the 
height range 40-70 km, and winds to lower altitudes. Modifications 
required in the 1965 data reduction program in order to achieve 
ROBIN measurements above 70 km have been described by Engler (ref. 8) 


227 



and Luers (ref. 2). Figure 1 indicates in greater detail, accord' 
ing to the author's count, the quantity of data obtained. 


This discussion of the data is from the standpoint of the user. 
In order to ascertain the usability of sphere data _ in such tasks 
as the construction of synoptic maps or the analysis of small sea es 
of motion, some review of the data accuracies is needed. It is 
stressed that odd features encountered in various soundings have 
been more a matter for explanation than alarm. A few soundings, 
however, have nearly defied explanation. The discussion of accur 
acies will therefore be followed by examples of unusual soundings 
which could be shown to be valid and some which remained suspect . 

The latter are rare, but they have required special probing by the 
user in order to judge their acceptability. 


THE VALIDITY OF ROBIN WINDS AND DENSITIES 


Our experience with ROBIN data indicates that with rare 
exceptions, the winds and densities are broadly representative of 
the ambient conditions. The rare exceptions might include condi 
tions of very large vertical motions (neglected in the drag equa 
tion for density) , or appreciable error in the drag coefficient, 
or unusual balloon behavior perhaps undetected by Engler s Lambda 
check (see below) (Jones and Peterson, ref. 3). With regard to 
small-scale variations of the wind and density, the fidelity of 
representation must depend greatly (as with any observational 
method) on the smoothing performed on the original data points. 


Criteria for Acceptability 


All soundings considered in this article meet the stated 
criteria for valid data (Engler, ref. 1). These are, as we have 
understood them: 


(1) Densities are within the stated accuracies, to be cited 
below, for that portion of a sounding satisfying the Lambda check. 
(Mathematical symbols are listed in the appendix.) Lambda is a 
measure of the vertical density gradient (A-p dP/dz) , e 
check assumes that reliable densities are obtained if Lambda, in 
practice approximated from the vertical acceleration data, falls 
within a defined neighborhood of the standard atmosphere value. 

Any unusual perturbation in the vertical sphere motion is assumed 
to be due probably to collapse or drastic change in shape o e 
balloon or perhaps to some unusual aerodynamic behavior. Wit 
non-spherical balloons, the C and cross-sectional area would not 
be known and the density coula not be ascertained. 


(2) When the Lambda check indicates collapse below 50 km 
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co f| sidered highly reliable. However, if the 

un^elLh^ Ck ' failS 50 - km ' the densit Y d ata are considered 

reliable, in accordance with the inference that the balloon was 
never fully inflated (.Lenhard and Kantor , ref. 9). aj - LOOn was 

. , , ^ Winds are considered broadly reliable at all times after 

Fot a a ioH- ha H aC ?? leratad t0 ValueS ° f 2 neater than -3 m sec^ 
For a non-rigid balloon, however, the response to the wind is not 

eJtimst^ asc ® rt * inab £ e from theory (ref. 9), and different error 
estimates apply (see below) . 

(4) Temperatures are considered broadly reliable at two 
scale-heights or roughly 15 km, below the starting altitude <^ee 
special section on temperature, below). 


Smoothing Interval 

f ROBIN data in the University of Dayton printout are given 

but what Y do e S? nd ab ° Ve 50 Km and f ° r ever y 2 seconds below 50 km, 
tiona? L?* ^ bese represent? Starting from 0.1-sec. radar posi- 

ai Sf the n\ 65 avera 9 ed to 9 et 0.5-sec. positions. Next, 

ff °£ the 5-sec. positions are fitted by one-degree polynomial 
least squares to give 15-sec. values for velocity (from the slSpe 

above G 50 U and ? 6 J e ^ ocities are recomputed at one-second intervals 

above 50 and 2-second intervals below 50 km, by dropping and addina 

data points at top and bottom; and then accelerations are determined 
by least squares fit to 7 of the velocities. This results in Ic- 
celerations valid for 22-sec. intervals above 50 km and 28-second 
£^ e £^ al ? be £ ow 50 km * For densities, which are proportional to 
review Jo era ?£° nS \ ! ame time . intervals apply. To make this 
solvina for^h^ 00 ^ 1 ^ posaible ' the pertinent equations for 
appendix f ° r ^ Wind ^ thermod y namic data have been stated in an 

interra^nfb^ veiocities are effectively determined over a time 
W? f 1/4-minute and the accelerations and densities over 
neariy 1/2-minute. For typical fall rates of the ROBIN one-meter 
balioon, a time interval of 1/4-minute corresponds to a descent of 
p n ^ m abo ^ e ^ km a nd less than a ki i ome ter below about 50 km. 

Engler (ref. 1) indicates that the horizontal distance of the 

reja?d2d as ref i lllatlons ? f peri ° d exceeding 1/4-minute, which he 
regarded as real His analysis of the effect of varying the number 

vfl show^f/S positions used for the basic smoothing inter- 
val shows that differences up to several m sec" 1 are possible in 

^Fia^pf^M 6 ° f thS °f cillations ' according to the data fit used 
(Fig. 2) More recently, Boer and Mahoney (ref. 10) have analyzed 
research series of ROBIN soundings for March 6, 1965 (White Sands) 
smoothing over a constant-height interval rather than a constant ' 
time interval. Various thicknesses were tried, from 100 to 500 
meters; these are substantially narrower layers than the effective 
smoothing interval at the higher altitudes of Engler 1 s data. For 
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these layers , Boer and Mahoney found that if the acceleration term, 
zx/(z-g), was included in the wind computations, the correlatio 
between wind profiles based on data from two radars (FPS 16) track 
ing the same balloon was greatly reduced. 

The reduction in one sounding (No. 933, 1500 MST) was so dra 
matic that it occurred to us to correlate the wind prof lies phased 
on the Engler data, which also include the acceleration or respo 
correction, but involve a smoothing interval of 22 sec For the 
actual fall rates of this sounding, this amounts to a layer thic 
ness ranging from 1.8 km to 3.4 km at 60 km. Figure s ows 
wind prof iles (u component), which are indeed in cl ° se . ” 

ment. The correlation between them, on removal of the basic trend, 
is 0.93^ in contrast to an extremely low correlation, 0.09, cal- 
culated by Boer and Mahoney . 

Thus, insofar as it rises above the radar noise level, Engler' s 
smoothing appears more realistic. At altitudes of 45 60 . 

Lettau (ref. 11) appears to have made effective use of the White 
Sands Engler-reduced data in tracking small-scale structure which 
he interprets as evidence of internal gravity waves. At high 
altitudes Engler 's smoothing interval, however, becomes quite gross 
(effectively, 4-5 km above 70 km), and consideration should be 
given to a shorter time interval which would exclude tracking error 
and yet permit the resolution of small-scale wave structure. o 
only would the time interval be critical, but the choice of poly 
nomial fitted to the high-altitude data, whether cubic or linear, 
for example (Luers , ref. 2), would also be important. In ee , l 
appears that various tradeoffs would be necessary to minimize the 
error in both wind and density, while suppressing radar error. 


Mention should be made of small-scale oscillations in the 
ROBIN density profiles. Unfortunately, there were very few valid 
thermodynamic data from the White Sands series. Figure 4 is a 
plot of the two profiles based on radar tracking of the same 
balloon (sounding no. 933). The curves are in excellent agreement. 
The indicated oscillations in both curves are of small amplitude, 
making a more detailed comparison difficult, but the oscillations 
are clearly in phase and with scarcely detectable divergence. Like 
the wind profiles, these data indicate that for the smoothing 
interval used by Engler, there is no apparent distortion from 
radar error. 


lcorrelatio n is meaningful in first decimal digit only, owing 

to subjectivity in determining the trend and choice of sampling 
frequency in the profiles. 
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Accuracy Estimates 


Finally, we need to have in mind the overall error estimates 
in the ROBIN data, as given by Engler (ref. 1). These are, for an 
FPS-16 or similar radar and the curve-fit (31-7) commonly used: 

TABLE I. ROBIN error estimates 9 






RMS 

Error 




Above 

60 

km 

50-60 

Below 

WIND 

Rigid balloon 

6 

knots 

2.5 

1 


Non-rigid 

10 



7 

5 

DENSITY 

Rigid balloon 

3.. 

5% 


3% 

3.5% 

PRESSURE 

Rigid balloon 

3% 



3% 

2.5.% 

TEMPERATURE 

Rigid balloon 

10% 



3% b 

4% b 


TEMPERATURES FROM SPHERE SOUNDINGS IN GENERAL 


The temperature error requires special discussion. As is 
well known, a relationship based on the hydrostatic equation is 
used for deriving the temperature when no thermodynamic data other 
than a density profile are available (eq. 3, Appendix) . An initial 
guess of the temperature T , is required at the starting altitude, 
i.e. the top altitude with°density data. The error in temperature 
is thus a function of two factors, (1) the departure of the initial 
temperature from the true temperature, and (2) the error in the 
density throughout the layer of integration. 

For an error-free density profile, it is evident that nearly 
ambient temperatures are not achieved until the ratio p /p be- 
comes negligibly small. This happens typically at an altitude 
roughly two scale heights, or about 16 km, below the starting 
altitude . The reduction of the ratio p /p with increasing 
height separation is shown in Table II. &or example, for a height 
separation of 2 scale heights p /p = 0.135, and the temperature 
error ranges from 0.5 to 3% depending on the error in the temper- 
ature guess at the starting altitude. The temperature error itself 


a 

Based on Table I (ref. 1) . 


b 

See discussion of temperature error in following section. 
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was obtained by evaluating the term, ^ AT, where A T is identi- 
fied with the potential error in T (eq. 3) for a given climatolog- 
ical regime . o 


If a certain level of temperature accuracy is desired, say 
2.5%, it is thus not generally possible to estimate the first 

good" data without some preconceived idea of the tem- 
perature variability at the station of interest. In low latitudes, 
we know from rocket grenade observations, for example, that if the 
temperature in the U.S. Standard Atmosphere Supplements, 1966 

12 ^ • ■‘■ s used for T n f the true temperature will not likely 
differ by more than 20°.° Thus, from Table III it can be seen that 
nearly ambient temperature would be achieved at a height separation 
o about 10 km. In high latitudes in winter, however, the tempera- 
ture variability is typically much greater, and it is doubtful that 
ambient temperatures will be sensed at separations less than 2 scale 
heights. The summer mesopause is known to have very cold tempera- 
ures associated with it, with the scale height possibly as low as 
4-5 km, so that near 80 km real temperatures might be sensed in 
relatively shbrt order. 


The other error source for temperature is error in the densi- 
^ es themselves and it should be noted that it is the error over 
the layer of integration that matters, not just the error at al- 
titude (eq. 3, Appendix). According to Engler (ref. 1) various 
density error profiles are possible. However, it may not always 
be possible to describe the height configuration of the density 
error in individual profiles, so that in some cases the total 
error m temperature may not be ascertainable. At the very least. 
Table III clearly illustrates that the temperatures provided in the 
first few kilometers below the first level of density data should 
not be construed as real temperatures. At times the reported 
temperatures may fortuitously come close to the real values, but 
we know of no way to readily distinguish these cases. The practice 
?* P ,! lshl ? g complete temperature profiles in the data books of 
the Meteorological Rocket Network thus seems questionable, at 
least without some qualification as to the validity of the data 
at the topmost levels. Jones and Peterson (ref. 3), however, 
consider publication justifiable on the grounds that although the 
absolute values of the temperature may be in error, "valid trends 
can often be seen." - 


For the height range in which reasonably accurate temperatures 
® . ®^P®cted , comparison with data obtained by other techniques 
is desirable. Since this topic was to be considered in depth by 
other speakers, only a few remarks will be made here. 

Jones and Peterson (ref. 3) have discussed the extent of 
agreement of data obtained with the aid of the University of 
Michigan 66-cm sphere, with grenade (layer-average) temperatures 
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and with HASP rocketsonde data in August 1965, Jt Wallops Island 
Initially, the sphere and grenade measurements differed appreciably, 
toe spiers tompe?atures being warmer. Improvements an the macro- 
Dhone array for the grenade measurements and in the sphere drag 
coefficients ?esulted in better agreement. In a combined sounding 
in which grenade and sphere measurements were separated at most by 
10 minutes, the temperatures were found to agree generally within 
a fei degrees. Oscillations in the sphere temperature curve were 
not presin? in either the grenade data nor in the HASP rocketsonde 
data at lower altitudes. Certain oscillations should,* of cours , 
not be expected in the grenade temperatures, m view of the layer 
averaging? their absence in the HASP data, however, suggests that 
™ oscillatory part of the sphere data may be erroneous. Jones 

and Peterson gave several possible explanations, name the°sphere 
radar track, a deflated sphere, peculiar aerodynamics of the spher , 
or the effect of large vertical motions in ^ , at ™° sphere . Th 
first two were not judged to be the cause in this case there were 
arguments against aerodynamic behavior as a cause; a ^d the' v 
cal motion effect 2 could not be evaluated owing to a lack of suit 

able information. 

In the case of temperatures from ROBIN soundings, S. Teweles 
(private communication to N. Engler, Jan. 7, 1964) pointed to an 
apparent discrepancy between ROBIN and grenade temperatures but 
subsequently Engler determined that the ROBIN data used for compari- 
son did not Y meet the criteria for acceptability (balloon collapse 
above 50 km) . Indeed, even now, there is no extensive set o 

robin and grenade data available, to ou ^ com-^ 

similar observing conditions, which would permit definitive com 
oarison The situation with respect to rocketsonde thermistor 
dato dSis no? seem much better, since in the altitude region where 
ROBIN temperatures should be most reliable, about 5 
ing the thermodynamic data commence at n.70 km) , the thermistor 
temperatures are subject to increasing error with height. The 
results of comparisons of ROBIN and Arcasonde temperatures at 
altitudes near 50 km will be presented m the next section. 


2 The subject of vertical motions was to be considered by 
another speaker? it is generally agreed that vertical motions ™ 
important only if they are of the order of m sec . A ^^ g 3 Q P ^ 
haps rare, motions of this magnitude may be possible_, above 30 km. 
in the stratosphere vertical motions of cm or mm sec 1 have been 
generally found (Miller, ref. 13), but during a stratospheric x 
warmina in 1966, upward motions as high as one-half meter sec 
were estimated (Quiroz, ref. 14). Above the mesopause motions near 
10 m sec" 1 have been reported. Thus the effect of vertical motions 
may be a difficult problem to evaluate until better statistical 
knowledge of the motions is obtained. 
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ILLUSTRATIVE DATA 


Cape Kennedy, Dec. 7, 1964 

A ROBIN sounding for Dec. 7, 1964 has been selected to illus- 
trate the broad agreement of the ROBIN data with atmospheric measure 
ments by other methods (Fig. 5). On this date the zonal wind 
indicated by the sphere increases at an extraordinary rate between 
30 and 55 km; the local vertical shear near an altitude of 33 km 
is about 20 m sec - -*- km - -*-, a value which may be considered statisti- 
cally rare. At 55 km the westerly wind component exceeds 100 m 
sec - -*- and is three times as great as the tropospheric wind maximum 
at the same station. So unusual a sounding might be questioned, 
especially in view of the westerly firing angle (to the east) at 
Cape Kennedy. 

Later the same day, however, a wind profile based on the 
tracking of a parachute indicated close agreement with the ROBIN 
profile (the observed differences are probably due to real atmo- 
spheric variation) . Figure 6 shows the synoptic situation on 
Dec. 9 (maps Were not analyzed for Dec. 7) ; a very strong jet is 
found at the 0.4 mb level (about 55 km over Florida), with strong 
winds observed as far west as Hawaii. At the 5-mb level (about 
35 km) , the winds are light over Florida and to the south, in 
agreement with the data for Dec. 7. 

The temperature and density profiles for this date are also 
of some interest. In view of our earlier discussion, realistic 
sphere temperatures would be expected at Kennedy some 10 km below 
the starting altitude and indeed, good agreement with rocketsonde 
thermistor temperatures can be seen just above 50 km. Near 55 km, 
comparison is precluded by the likelihood of increasing error 
with height in the thermistor measurements (Quiroz, ref. 15). The 
reality of the temperature difference at 42-44 km would be difficult 
to ascertain. It seems reasonable to conclude that reliable 
temperatures are indicated by both methods of observation at least 
in the height range 44 to 52 km. A statistical comparison of 
sphere and thermistor temperatures will be described below. 

The ROBIN densities are 7-10% lower than the Arcasonde 
densities, for which the observation time is in mid day. Part of 
this difference may be due to diurnal variability. The possibility 
of a systematic bias will be explored below. 

Ascension Island Temperature Comparison, March-June 1964-65 

The period March-June 1964-65 was chosen for the comparison 
of temperatures from Arcasonde thermistors (1965) and from ROBIN 
sphere soundings (mainly 1964). The comparison was limited to 
three altitudes, 46, 50, and 54 km, where data from both sources 
should be considered reasonably reliable. The thermistor data have 
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not been corrected for possible aerodynamic, radiational, and con- 
duction heating errors , which have been estimated at about 2 de- 
grees (total) at 46 km, increasing to about 5 degrees at 54 km 
(Drews, ref. 16). The number of ROBIN soundings (18) is small, but 
nevertheless represents one of the densest clusters of such sound- 
ings available, save for the experimental series at Eglin AFB in 
1960-62. 

The observed temperatures are plotted in Figure 7 and relevant 
statistics are given in Table III. Inspection of this table shows 
that: 


(1) Average temperatures from the two sources agree within 
about three degrees if the thermistor data are un _ 
corrected . 

(2) Sphere temperatures are warmer by a few degrees if 
compared with corrected thermistor data. 

(3) The dispersion of the sphere temperatures, as given by 
the values of standard deviation, is greater than for 
the thermistor temperatures. (No attempt was made to 
smooth the oscillations present in the sphere data; 
smoothing would have brought the standard deviations 
into closer agreement.) 

While the sample is probably too small to give stable statis- 
tics, these data indicate that over a definable height range the 
sphere average temperatures are at least realistic. Further com- 
parison in a regime of greater variability (middle or high latitudes 
in winter) is desirable. Moreover, the influence of oscillations 
of large amplitude needs to be examined further. 

An interesting feature in Figure 7 is the indication of a 
diurnal temperature increase from 04-05 GMT to 16-18 GMT in three 
pairs of observations. At 46 km, the two thermistor pairs on 
May 23 and 26 indicate a diurnal range of about 10 °C; this range 
is also indicated by the pair of sphere observations on April 8, 
which are not subject to any direct radiational error. At 50 km 
a similar behavior is observed? part of the large temperature in- 
crease in the sphere pair, however, may be due to non-diurnal effects 

Enigma in Ascension Island Density, August 1964. 

Figure 8 depicts the observed densities at Ascension Island 
in 1964 at two altitudes, 46 and 60 km, based on ROBIN sphere 
soundings, together with comparative data obtained by other methods. 
The lower altitude was chosen because at this height the sample of 
thermodynamic data from descending spheres was still appreciable 
and would permit comparison with values derived from thermistor 
measurements; as previously indicated, the error in the latter 
has been considered small in the upper stratosphere. 
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Figure 8 shows a number of interesting features, of which 
the most striking is the unusual and highly improbable behavior in 
mid-August. A density increase by 35% at 46 km and by 43% at 60 km 
is indicated by the sphere data in a period of 31 hours. These 
values greatly exceed the maximum density change in 24 or 48 hours 
previously indicated from rocketsonde data, namely a change of 19% 
occurring at a high— latitude station (Quiroz, Lambert, and Dutton, 
ref. 17). Yet by the criteria for acceptability stated earlier, 
the soundings in August must be considered reliable soundings . 

It is noteworthy that the extraordinarily high density on August 17 
was again observed two days later, and the low value of August 16 
was similar to the value observed earlier on August 12. Corre- 
spondence with the Air Force office responsible for the observation- 
al program at Ascension Island did not reveal any irregularity in 
data reduction. Subsequently, Engler (ref. 1) applied his time 
of-f all test for suspect balloons and judged the soundings to be 
valid. Thus, while the values observed on August 17 and 19 were 
too extreme to inspire credibility, there seemed to be no way of 
showing that the data were incorrect. In preparing this review, 
it occurred to us to examine, insofar as possible, the internal 
consistency of the thermodynamic data with the observed winds. 

The temperature change from August 16 to 17 amounted to only a few 
degrees at 46 km. Thus a large pressure increase was associated 
with the increase in density. Fortunately, pressure data were 
also available for another rocket station, Antigua (17°N) , on 
August 17. A geos trophic computation, assuming a linear pressure 
change between the two stations^, indicates that if the pressures 
and densities at Ascension are valid, the zonal wind at Ascension 
should exceed 400 m sec - ^-. The wind obtained from radar tracking ^ 
of the ROBIN sphere was at most 25 m sec -1 in the vicinity of 46 km. 
We therefore conclude, on the basis of purely physical reasoning, 
that the soundings of August 17 and 19 are invalid. Rocket grenade 
observations were also taken on August 16 and 17 , though the 
results were not available until much later (Smith et al., ref. 18) . 
These data, entered on Figure 8, are completely at variance with 
the sphere results. 

This case is but one, although possibly the most dramatic, 
of several extraordinary sphere soundings encountered by the author, 
and it emphasizes the need for careful scrutiny of all data by the 
user. Indeed, all rocket soundings by whatever method require 
careful review, since the many aspects of data reduction and 
transcription may increase the possibility for error. 


A non-linear pressure distribution would require an even 
greater wind than that computed, at some point between the two sta- 
tions . 
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Densities from Sphere Soundings Compared with Other 

Measurements 


Another feature of interest in Figure 8 is an apparent dis- 
crepancy in densities based on sphere soundings and on thermis- 
tor measurements. Ascension Island offers a useful opportunity for 
comparison because at the end of September 1964 the observational 
program reverted from an exclusively sphere schedule to a predomi- 
nantly Arcasonde schedule. Also available for comparison in 1964 
are a few grenade and Pitot-static tube measurements at Ascension 
and several University of Michigan sphere soundings at Kwajalein 
(9°N, 168°W) . From other years, tropical grenade soundings at 
Natal, Brazil (6°S, 35°W) (1966-67) and Guam (14°N, 145°E) (Nov- 

ember 1958) are also entered. 


Two points are readily apparent: 

(1) The sphere densities are, with rare exceptions, lower than 
the mean based on thermistor densities; e.g. in September-December , 
1964, the mean sphere density is nearly 10% less than the mean 
based on thermistor measurements. 


(2) A greater dispersion is indicated by the sphere measure- 
ments . 


At the upper altitude, the difference between sphere and ther- 
mistor measurements is partly due to the error in the latter, which 
increases strongly above 55 km, and no attempt has been made to 
enter individual thermistor values for 60 km. 

It is therefore meaningful to concentrate on the data for 46 km, 
where the error in either set of data should be minimal. Various 
possible explanations of the observed difference merit considerations, 
such as unsuspected temperature error in the rocketsonde measure- 
ments, error in the drag coefficient used for the sphere reduction, 
etc . 


The effect of a thermistor temperature error may be evalua 
with the aid of the integrated hydrostatic equation in the form 



4 



( 1 ) 


In practice, pressures and densities based on rocketsonde 
temperature are obtained through use of the approximation 
p = p o exp (-gAz/R T) , but as has been shown by Ballard (ref. 19) , 
the departures from results based on the more exact Eq. (1) are 
negligibly small. 
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Introducing an error, AT, in temperature, held constant over the 
layer of integration, the density at the upper altitude becomes 



and the ratio of (2) to (1) is 



For a temperature structure approximating the conditions at 
Ascension Island 5 and for a hypothetical temperature bias, AT, in 
the thermistor soundings, sustained over the 20-km layer from 26 
to 46 km, the solution of Equation (3) yields a density error of 
approximately +3% if AT = +3°, increasing to +8% if AT = + 8°C. 

(Over a 10 -km layer from 36 to 46 km, a temperature error of 16° 
would be required to explain a discrepancy of 8% in the density at 
46 km.) The rocketsonde temperatures are generally believed to be 
quite accurate below about 50 km, although some disagreement with 
radiosonde temperatures near 30 km has yet to be explained (Quiroz, 
1969) . Since an unreasonably large temperature error is required 
to explain the density difference in Figure 8, it appears that we 
must look to some other error source, or more likely, a combination 
of sources. 

According to Engler (ref. 1) uncertainty in the drag coefficient 
for ROBIN spheres is less than about 2% below 50 km, but improved 
knowledge of this factor is needed. Luers(ref. 2) and Peterson 
(ref. 6) have pointed to inconsistencies in the available drag tables 
and have re-emphasized the need for improved data. An interesting 
series of measurements with hypersonic rigid spheres (Kwajalein, 
1965-1968) has been obtained under conditions for which a high 
degree of confidence can be placed in the drag data used, according 
to Salah (ref. 4, 21). Comparative ROBIN and hypersonic sphere 
measurements, if feasible, might shed light on the drag data used 
for the inflatable, subsonic spheres. 

It is interesting to note in Figure 8 that measurements by the 
grenade method tend to lie between the sphere and thermistor values. 
Almost without exception, the grenade densities are lower than the 
thermistor values. It is therefore our belief that the discrepancy 
between the sphere and Arcasonde data must be due to error in both 
methods of measurement. 

5Supplemental Atmosphere data for 15°N were used. 
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Diurnal series, Eglin AFB , May 9-10, 1961 

series n ?“ e LvTlJ e Sp !?? re data is found ln the ai«rnal 

indicates a Xn = 7)°' ' 96 (Flga re 9 >- The 2 °°° GMT observation 
a d f nsit Y increase of nearly 30% from 1600 GMT. This 

Inrt f T s ° undin ^ s at 1015 and 1115 were not used by Cole 

and Kantor (ref. 20) in their harmonic analysis for data at 60 km 
Summing the first two harmonics. Cole deduced a daily range of 12% 

par^ulJ^-d J h V nC ^i° n ° f the three unused sS s in 

particuiar the data for 2000 GMT, would probably have had a strona 

Save not 8 ™? ^ r “ ults - T ° the author's knowledge, these data * 

n °t yet been shown to be invalid. Lenhard's (ref. 22) analysis 
of the wind data for this observational series indicated a weSkeni™ 
' ? aster jY flow in Parallel with the afternoon rise in lensUv 

?eauf?ed"?o a t 1Urnal ? dveotive <*««*, but careful study would be 
required to show a relationship. 

USES OF SPHERE DATA AND FINAL REMARKS 

var i‘ ptv^nf fr ° m falling spheres have already proved useful in a 
They have already been used S) P ° tential has not been exploited. 

(1) the construction of high-level synoptic maps. 

v . , (2) Preliminary determinations of diurnal variability of 
density and wind: m the lower mesosphere from ROBIN data in the 

?ef 2^ he Jnd a ?n "V es ° s P h -e.fro m Australian spheres (Ro^e et aU 
? d ^ he < 3 uasi-ls °Pycnic layer at about 90 km fFoirT 
Michigan spheres (Jones and Peterson, ref. 3) . 

M V, exploratory studies of small-scale variability (Newell 

Mahoney, and Lenhard , ref. 24, Lettau, ref. 11; MaJo^ey and loer 
ref. 25; Cole and Kantor ref. 26 and others). ' 

90-120 (4 km m ?n e J S Sf t l 6 ? e !?f ity structure in the important region 
yo 120 km (U.S . Standard Atmosphere Supplements, 1966, ref. 12). 

Salmela^n^Q^ 010 ^ 0 ^ da I: a P rocessin 9 (e.g., Brockman, ref. 27, 
baimeia and Sissenwme, ref. 28). 

With regard to (1) , the synoptic analysis program of the 
Upper A lr Branch, National Meteorological Center, has resulted in 
a contrnurng series of weekly constant-pressure iharts a? levels 

“ f a )const ant ^ l 42 < a nd 55 km (5 ' 2 ' and °- 4 "*>, beginning 
19 ui (Constant- level density charts have also been produced 

though less frequently.) Because of the preponderance of rocket- 

sonde data (Quiroz, ref. 15), the utilization rltltlf sphere 

for iitut^nd n t 6 analysis technique requires wind and^emperature 
for input, and not surprisingly the temperatures from ROBIN spheres 
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could be used only infrequently. This was due to (J'^he typically 
mrill Vic. i rrh t- ranae of valid temperatures and (2) untenable depart 
” f 5 e otherwise smooth temperature fields depicted an the 
analvses These departures were at times associated with oscilla 
torv Y features in the sphere profiles and would have been minimized 
if smooth profiles had been used. Nevertheless, the sphere data 
ifveaf times provided missing links in a rocket network of sparse 
coverage , and ?hey have a strong potential utility in future work, 
particularly if the height range of valid data is extended. 

One of the most promising uses of the data oo^iT^Only a 
defining the complex density structure from 90 . 

handful of soundings have been used as a basis for en g 

structural models of the thermosphere, yet as Figure 1 shows there 
taan a larae increase in the number of soundings obtained in 
recent years. In addition to Priding a better grasp of boundary 
conditions for thermospheric models (Thomas, ref. 29), P 
knowledge will permit more definitive investigation of the inf lu 
ence o^variable solar activity at these altitudes (Lmdblad, 
ref. 30; Ellyett, ref. 31). 

We have sought to indicate, through a few examples, that the 
sphere density and wind data are, with rare exceptions reliable 
nLr definable ranges of altitude. Temperatures should be reliable 
at G two scale heights below the first altitude of density data, but 
the oscillations in temperature, which either are of ^^ater amp 
tude than those encountered m rocketsonde prof lies or 
times not present in the latter, need further stuy. 
ommended that this problem be given special attention, since tne 
full utility of the temperature data cannot be achieved unti 
reality of the oscillations is ascertained. Other problems for 

investigation are the apparent systematic dif f^accuScy 6 of 
based on sphere versus thermistor soundings, the accuracy of tne 
avSaUe d» dlta, and the possibility J«ge -rtical muttons 
which miaht at times affect the sphere results. With regard to 
hicrh-altitude data from spheres with fast descent rates, consider 
ation should be given to determining an optimum smoothing interval 
SmcE Tt tie same time suppresses radar error and preserves small- 
scale atmospheric structure. Finally, it is recommended that 

comparative experiments in the future be condu ^ ed a ?or the inam- 
winter in high latitudes, under conditions which favor the unam 
biguous separation of observational error from true variability. 
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appendix 


PHYSICAL EQUATIONS FOR REDUCTION OF SPHERE DATA 


Atmospheric 

Variable 


WIND 


Complete 

Equation 


Simplified 

Equation 


(z - w) 


u = x - 


( s 


+ C - B 1 
x xl 


(la) 


U = X - 


z + C - B 
z z 


z x 


Z - g 


(lb) 


u, horizontal west-east component 
w, vertical motion (atmospheric) 
g, acceleration of gravity 
C , B Coriolis, buoyancy forces 

Equation for v, the horizontal south-north component, 
is similar. The total wind, V = (u + v ) / . For 
generality, the Coriolis and buoyancy forces m (la) 
are retained in the ROBIN program. C is judged 
significant above 90 km (ref. 1). The buoyancy force 
B = (v) B pg, may be large below about 30 km. In the 
simplified equation (lb) , vertical motion and 
horizontal Coriolis and buoyancy forces are neglected, 
This eq. is used after collapse of ROBIN balloon. 

Under certain conditions, the acceleration or 
"response" term, z x/(z - g) , might be neglected if 
u = x, but Engler shows that error due to neglect of 
this term may be very large at high altitudes and 
does not recommend the use of the approximation. 


DENSITY In general, F D = l/PpVgCpA = ma D 

2maj) 


i . e . , 


P = 


V | C D A 


For ROBIN 


P = 


a - C U 

CjjAVU - w) + (v)gg z 


p = 


-2m(z ~ g) 


C d AVz 


(2a) 


(2b) 
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APPENDIX (Continued) 


In (2b) , Coriolis and buoyancy forces and vertical 
motion are neglected. This is form of eq. commonly 
used . 


TEMPERATURE 


In general , 
Integrated 

T 


Sp/3z = - pg; p = pRT . 
hydrostatic equation, 

1 r z ° . p„ 

' # J z Pg + T? T o 



Pg Az + 



T 


o 


(3a) 


(3b) 


z q is starting altitude. Eq. (3b) is used. 

Mathematical symbols not defined above: 
x, x dx/dt, d 2 x/dt 2 , etc. 

A sphere cross-sectional area 

a p drag acceleration 

C d drag coefficient 

Y temperature lapse rate, dT/dz 

H pressure scale height (H = RT/g) 

m sphere mass 

p pressure 

R gas constant for dry air 

P air density 

Srj, standard deviation of temperature 

T , T temperature, mean temperature 

V B sphere velocity 

(v) sphere volume 
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altitude km 


64 



Wind speed m sec 


Figure 2. -Variation in amplitude of indicated 
due to choice of data fit. (Adapted from 
ret. l, which shows additional curves for 


wind oscillations 
Figure 9, Engler, 
other fits.) 


ALTITUDE KM 


ALTITUDE KM 



Figure 3 . — Plot of wind profiles 
radars tracking same sphere 
Engler smoothing interval. 


determined from data of two 
, indicating high correlation for 


ALTITUDE KM 




DENSITY g m 3 x 10 



Figure 4. -Density data for same sounding as in fig. 3, showing 

excellent agreement in values based on data from two radars. 




launching 6 hours later. 
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Figure 7 . -Comparison of temperatures from sphere soundings, 1964 
65 , and thermistor soundings, 1965 , at Ascension Island. 

See statistics in table III. 
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Figure 9. -Densities indicated in May 1961 diurnal series of 

sphere soundings. Circled values were apparently not used 
in harmonic analysis of 60 km data by Cole and Kantor 
(ref. 20) . 
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Departure from mean 



